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Abstract. - Mixing of two fractions of a granular material in a slowly 
rotating two-dimensional drum is considered. The rotation is around 
the axis of the upright drum. The drum is filled partially, and mixing 
occurs only at a free surface of the material. We propose a simple 
theory of the mixing process which describes a real experiment sur- 
prisingly well. The dependence of the mixing time on the drum filling 
is calculated. The mixing time is infinite in the case of the half-filled 
drum. We describe singular behaviour of the mixing near this critical 
point. 

A beautiful experiment on so called avalanche mixing of granular solids 
in a slowly rotated two-dimensional drum was presented in the paper of Met- 
calfe et al [|l|]. The patterns observed were so striking that one of the figures 
made the cover of the March (1995) issue of Nature. In the present com- 
munication, we shall propose a simple theory which provides a surprisingly 
good description of this experiment. We shall demonstrate that kinetics of 
avalanche mixing can be understood using a geometrical approach without 
appealing to ideas of self-organized criticality f2|. (One should note that the 
flow of granular solids in rotated drums is usually used to study self-organized 
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criticality §, |, § |, % §, g [TO|, [IT], |IJ, etc.) In fact, we shall find the 
main quantity appearing in the problem — the mixing time — for all values 
of the drum filling. 

Two types of small granules - - black and white — are poured into a 
two-dimensional upright drum. Except of their colours, the granules of the 
different fractions are indistinguishable. The drum is filled partially, so there 
is a free space at the top. The specific initial distribution of pure fractions will 
be insignificant for our aims. (In the experiment at the first moment, the 
black fraction is above the white fraction — see fig. 1, a.) The drum starts to 
rotate adiabatically slowly around its axis in a counter-clockwise direction, 
for definiteness. Therefore, the complete angle of rotation plays the role of 
time. The answer will not depend on the radius of the drum, which may be 
assumed equal to 1. It is convenient to characterize the relative volume of the 
unfilled space by the angle 6 (see fig. 1, b). We describe macroscopically the 
state of each point of the material by the quantity p — the concentration of 
the black fraction at a given location (for a pure black fraction we have p = 1 
and for a pure white fraction p = 0). Thus, we shall study the evolution of 
the space distribution p(x, t). 

Let us assume that the granules cannot intersperse and slip mutually as 
long as they are in the bulk of the material. Only the granules at the free 
surface layer flow and undergoing mixing through the action of successive 
avalanches. Such mixing is called avalanche mixing |]]. In a real experiment, 
the free surface is flat, and when the drum rotates slowly, the surface is tilted 
at the angle of repose (strictly speaking, the angle fluctuates slightly with 
time U, but these variations are small, and we shall neglect them). The 
angle of repose does not appear in our answers, so we set it equal to zero in 
fig. 1, b. 

The described simple configuration demonstrates clearly the essence of 
avalanche mixing. One can see immediately that if the drum is more than 
half full, the central part of the material — it is surrounded by a dashed line 
in fig. 1 — rotates with the drum without any mixing, and mixing occurs 
only inside the ring cos 9 < r < 1. If the drum is less than half full, the whole 
of the material will mix. 

Avalanche mixing proceeds in a fastest possible way, if the material is 
completely mixed after interspersion from the right half of the free surface 
to the left half (see fig. 1, a). For given filling of the drum, the mixing 
time for any material undergoing avalanche mixing cannot be shorter than 
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the mixing time obtained under this assumption. In this case, there may 
be inhomogeneous distribution of the fractions at the right half of the free 
surface, but the material becomes to be mixed homogeneously at the left 
half. After some time 5t < 2tt the material from the left half of the free 
surface turns as a unit to a definite tangent like that denoted as CD in fig. 
1, b. Therefore, at sufficiently long time t > 2{n — 9), the concentration p is 
constant at all points of each separate tangent to the dashed-line circle (see, 
e.g., the tangent CD in fig.l, b). On the other hand, the values of p may be 
different for the different tangents. Now it is possible to describe exhaustively 
the state of our system at the time t > 2(ti — 9) by the quantity p(f, t), where 
(p is the angle of the corresponding radius- vector with the normal to the free 
surface. Thus, in fact, the problem becomes to be a pure geometrical one, 
and there is no need to take into account all subtleties associated with the 
structure and shape of the granules. We shall consider kinetics of avalanche 
mixing in this fastest possible (for avalanche mixing) regime and compare 
the results with the experimental data |]J. As it will be shown, the regime 
under consideration is very close to a real one. 

Let us derive equations for the mixing dynamics. First we consider the 
case in which the drum is more than half full (8 < vr/2) and t > 2(n — 9). 
(Avalanche mixing at shorter times is discussed in our papers []14j| .) Let us 
turn the drum by a small angle dt at the moment t. Then the following 
amount of the black fraction will flow from right to left (from sector A to 
sector B in fig. I, a): 

Sin 2 9 rsinO 

dt——p{0,t)=dt drrp right (r,t) = (1) 
2 Jo 

, f 2S ,Jr(0 _ , _ , , cos 2 9 f 2S , A dtan 2 (C/2) tn . 
= Jo d ^^ r ^P^hMC),t) = dt^-j o d( J^ p ( 2 7r-C,t). 

The first integration Jq 1110 dr is from the free surface centre to the drum wall. 
r(C) is the coordinate of the intersection of the tangent fixed by the angle 
(f = 2it — ( with the free surface (see the segment ED in fig. 1, b). Note 
that r(() = cos6'tan(C/2), so r(29) = sin#. Thus, the concentration of 
the black fraction at the different points of the right half of the free surface 
equals p r ight(r((),t) = p(C,^)- One can see that p((p,t) = p(0,t — <p) when 
t > > 0, because granules cannot mix inside the material. Thus, from eq. 
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flip, we obtain the following equation for p(0, t): 



p(0, t) = cot 2 9 ^ d(p(0, t-2-n + Q sin(C/2) . (2) 
Jo cos^^/z) 

The case in which the drum is less than half full may be considered in the 
same way. After introducing the angle •& = 7r — 9 we obtain the equation: 

p(0,t) = cot 2 •& d(p(0, t - C) ■ (3) 

Thus, a complicate dynamical problem is reduced to studying of dynamics 
of the linear zero-dimensional system! Note that eqs. @ and (H) are not 
interconverted after the formal change of variables 9 = it — d because of 
the non-analyticity at the half-filling point 9 = d = n/2. Indeed, if the 
drum is exactly half full, it follows immediately from eq. (|2p or (|3p that 
p(0, t) = p(0, t — ir). Therefore the material will never be mixed to a uniform 
state, and the mixing time turns to be infinite. 

Let us study asymptotic behaviour of the solutions of eqs. (0) and (|3|). It 
is sufficiently to substitute p(0, t) = e zt into them. If the drum is more than 
half full, one gets from eq. (0): 

e 2 ^ = 2cot 2 # f e d^^f . (4) 
Jo cos 3 £ v ' 

If the drum is less than half full, it follows from eq. (|3p that 

l = 2cot 2 tf [\ e - 2 *-^- . (5) 
Jo cos 3 £ v ; 

We have derived two transcendental equations for z. To obtain the char- 
acteristic time of mixing one should find their nearest to zero roots. These 
roots are expressed by the following relation: 

z = -t~ x ± i2n/T , (6) 

where r and T are the mixing time and the oscillation period, correspond- 
ingly. Therefore, at long times, the concentration of the black fraction at the 
left half of the free surface evolves in the following way: 

p(0, t) = poo + Co exp{-t/r} cos {2vrf/T + <p } • (7) 
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Here poo is the concentration of the black fraction in the homogeneously 
mixed material, poo, Co, and <p are constants depending on the initial con- 
ditions, that is on the amount and distribution of the pure fractions at the 
initial moment. One can see that r and T depend only on one parameter, 
that is on the drum filling. Let us find the roots of eqs. (f|) and (||). 

Though the solutions of the transcendental equations are not expressed 
analytically, they may be found with any desired precision. The integrals in 
eqs. ([|) and (|5|) may be expressed in terms of special functions. Since these 
functions are rather exotic, we shall calculate directly the integrals in eqs. 
(Q) and (|5|). For example, one may integrate them by parts and then iterate 
the following relations: 




-2zcot 2 #^ ^e 22? tan 2 ^| (8) 

if the drum is more than half full (8 < 7r/2) and 

-4 log 1 1 + 2z cot 2 $e 2z <> die~ 2z ^ tan 2 f 1 
2v [ Jo J 



z 

2d 



(9) 



if the drum is less than half full < vr/2). 

The iterations converge very fast. Thus, using eq. (|^), we obtain the 
dependence of the inverse mixing time r -1 on the value of the drum filling 
(see fig. 2), that is our main result. We express r -1 through v = $ — sin 2$/2, 
i.e. the fractional volume of the drum which is occupied by the material. For 
comparison, in fig. 2, we also show the experimental points taken from the 
paper f|. 

The obtained r(v) is the shortest possible time of avalanche mixing at 
the given filling v. Thus, the curve t~ 1 (v) for any granular material must be 
under the extreme curve in fig. 2. One sees that though the experimental 
points from the paper |ij are near the extreme dependence, in fact, they are 
under it. We find that the experiment is very close to the regime of the 
fastest avalanche mixing. Therefore, it is described by our theory, which does 
not include any parameters of the mixing granular material. 

One may introduce the mixing rate v/r JTJ. In our theory, its maximum 
value equals max(f/r) = 0.0683. . . at v = 0.177. . . These values are also 
close to the experimental ones: max(v/r) = 0.056±0.006 at v — 0.17±0.015 
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The analytical results for r 1 and T for the limit cases of low and high 
filling are known already |f[4j| . For low filling t? -C 1 we have obtained r _1 = 
1.392/21? and T 4^^/7.553 and for the almost full drum (9 < 1) - 
r^ 1 ±^# 2 and T = 2tt — 48/3. 

Let us discuss the behaviour of r and T near the half filling point. If the 
drum is exactly half full, the concentration of the black fraction is a periodic 
function with the period n, i.e. the half period of the rotation of the drum, 
and the material will never be mixed. Thus, let us write the root sought in 
the form z = 2i + s where s is a small addition. First we consider the case 
in which the drum is more than half full. Substituting this relation into eq. 
(H) we obtain 

e 2 ™ = 2 cot 2 9 J 9 d£e 2sli {a sin 2£ + tan 3 f - 7 tan f - 4i(cos 2£ + tan 2 f - 1)} . 

(10) 

Now we introduce a small deviation 8 from the half filling: 6 = tt/2 — 8, 
and then one may expand e~ 2s ^ in small s£. Using the integration by parts 
it is easy to select the singular at 8 — > terms. The rest (finite at 8 — > 0) 
integrals are also evaluated without trouble. Then one should compare the 
coefficients of the terms having equal powers 8. 
After the straightforward calculations one gets 

z = 2l - -i - pog(l/|*|) - 5/2] + . . . (11) 

7T 7T 

Taking into account eq. (§), we find from eq. (|TTD that at sufficiently long 
times (t ^> t) mixing evolves according to the law (0) in which the charac- 
teristic mixing time r and the oscillation period T are expressed in terms of 
8 = n/2-9: 

_i 165 2 

T = 

7T 

and 

T = vr + 4<5. (13) 

r and T may also be expressed in terms of v — 1/2 = 28 /it. 

Similar calculations may be done for less than half filling. The results 
are the same as (|i~2"|) and ([U|) (now 8 = $ — n/2 < 0). Note that such 
a coincidence is unusual for critical phenomena Jy|. Analytical expression 
( P^[) for r(f) is available only in the narrow region \v — 1/2 1 < 0.02. 



log( 



(12) 
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In our paper we neglect fluctuations of the inclination angle of the free 
surface. In fact, this angle fluctuates in the interval between the angle of 
repose and the angle of marginal stability. If we shall take into account these 
fluctuations, the results will change slightly. For instance, in the limit of low 
filling v — > one obtains nonzero r. 

In summary, as we have shown above, in the experiment Jl|, mixing pro- 
ceeds in a nearly fastest possible (for avalanche mixing) way. In this case, the 
description of avalanche mixing becomes to be very simple — the problem 
is reduced to studying of the dynamics of a linear zero-dimensional system. 
In fact, we have solved it analytically. One see that the system under con- 
sideration demonstrates unusual critical behaviour near the half filling point 
at which the mixing time is infinite. Note finally, we have described quan- 
titatively the beautiful experiment without using any characteristics of the 
granules. 

The author thanks Yu. A. Firsov, A.V.Goltsev, S.A.Ktitorov, E.K.Kudinov, 
A.M.Monakhov, and B.N.Shalaev for many helpful discussions. This work 
was supported in part by the RFBR grant. 
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FIGURE CAPTIONS 



Fig. 1. - a) The avalanche mixing scheme. For an infinitely small turn 
of the drum, the granules of different fractions flow from sector A to sector 
B, undergoing mixing. The free surface of the material is tilted at the angle 
of repose all the time. The pure fractions are shown by the black and white 
colours. Regions of mixed material are denoted by the grey colour. The 
degree of mixing is not shown, b) The arrangement of the material for the 
case in which the drum is more than half full. The free surface is shown to 
be horizontal, since our results do not depend on the angle of repose. The 
angle 9 characterizes the relative volume of the free space. The position of 
the tangent CD is fixed by the angle (p. For the fastest avalanche mixing, 
the concentration of the black fraction at all points of such a tangent equals 
p(ip, t) at sufficiently long times. The tangent crosses the free surface if 
( < 29, where ( = 2n - (p. 

Fig. 2. - Inverse characteristic time of mixing as a function of the relative 
filled volume of the drum. For comparison, the experimental points from the 
paper || are shown. 
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Fig. 1 




Fig. 2 



